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Abstract 

We analyze the factorization process for lattice maps, searching for integrable cases. 
The maps were assumed to be at most quadratic in the dependent variables, and we 
required minimal factorization (one linear factor) after 2 steps of iteration. The results 
were then classified using algebraic entropy. Some new models with polynomial growth 
(strongly associated with integrability) were found. One of them is a nonsymmetric 
generalization of the homogeneous quadratic maps associated with KdV (modified and 
Schwarzian), for this new model we have also verified the "consistency around a cube". 



1 Introduction 



Although many interesting results have been obtained for discrete integrable systems and 
many properties have been clarified, we are still to great extent in the "taxonomic" stage of 
development. We do not have any clear classification and probably we have only discovered 
a small sample of integrable difference equations, the tip of the iceberg. 

As with differential equations, there is no universal definition of integrability for discrete 
systems, and consensus about integrability can exists only within certain subclasses of equa- 
tions. In this situation "integrability predictors" are very useful. These are algorithmic 
methods based on a somewhat weaker definition of integrability, which nevertheless seem to 
be associated with integrable systems. For example, for a Hamiltonian systems of 2N degrees 
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of freedom "Liouville integrability" means the existence of N independent and sufficiently 
regular functions that commute w.r.t. the Poisson bracket. Proof of integrability would then 
require the construction of the said commuting quantities, which is not algorithmic. In this 
case a good algorithmic integrability predictor is the Painleve test. 

1.1 Integrability test and concepts in ID 

For difference equations there is also a need for integrability predictors. One such method 
is "singularity confinement" (SC) [1], which has been advocated as the discrete analogue of 
the Painleve test. The idea of this test is to check what happens at a possible singularity 
of the evolution. Something special can only happen, if the value of the dependent variable 
becomes infinite. A singularity is then defined as a point where the next step cannot be 
determined (for example due to expressions like oo — oo). One then studies what happens 
near this point and the test condition is the following: If the dynamics leads to (or near) a 
singularity then after a few steps one should be able to get out of it, and this should take 
place without essential loss of information. 

This principle has been used successfully in deriving new integrable difference equations, 
especially in finding discrete analogies of Painleve equations in [2] and in numerous subse- 
quent papers by many authors, especially trough the so-called "de-autonomizing" procedure. 
Despite its success, it turns out that passing the singularity confinement test is not sufficient 
for regularity, counterexamples are given in [3]. 

It turns out that singularity confinement is strongly associated with reduced growth of com- 
plexity. When one iterates a rational map the expression becomes more and more complex 
as a rational expression of the initial values. If we measure this complexity by the degree of 
the numerator or denominator then generically the degree grows exponentially. However, the 
growth can be reduced if some common factors can be canceled. Some amount of cancella- 
tion always happens when singularity is confined 0], this is why SC is a such an efficient test. 
However, as far as integrability is concerned, it is the precise amount of cancellation that is 
crucial. This kind of complexity analysis, alias algebraic entropy calculation has turned out 
to be an unmatched integrability test for maps[T2| [3]. The conjecture about growth after 
cancellations and integrability is as follows: 

• growth is linear in n =>■ equation is linearizable. 

• growth is polynomial in n =>- equation is integrable. 

• growth is exponential in n =>- equation is chaotic. 

1.2 Integrability tests and concepts in 2D 

In principle the idea of singularity confinement can be applied to lattice equations as well, 
this was already discussed in paper pQ where the test was introduced. More recently this idea 
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was applied in [3] where an "ultra-local singularity confinement" was proposed (see also [0]). 
Nevertheless, it turns out that SC is rarely used in the study of 2D maps, perhaps due to the 
possibility of many different singularity confinement patterns depending on arrangement of 
initial values. In contrast the growth of complexity analysis can be easily applied to lattice 
equations as well (see (T5J EE]), and we will use it here. 

Perhaps the strongest form of integrability is "Consistency Around a Cube" as this kind of 
consistency immediately produces a Lax pair [8]. The idea here is that one should be able 
to consistently extend a two-dimensional map into three dimensions. If the 2D map is defined 
on an elementary square of the 2D lattice, then one constructs a cube with a suitably modified 
maps on all sides of the cube (the modification deals with associating different parameters to 
different coordinate directions). This creates a potential consistency problem in the evolution 
as follows: suppose we are given the values at the corners xooo, ^ioo> #oio> £001 (initial values), 
then the values at xhq, £101, ^011 ar e uniquely determined by using the proper maps, but the 
value at xm can be computed in 3 different ways. The consistency test is that these three 
ways should all yield the same value. This is a kind of Bianchi identity. It has been used as 
a method to find and classify integrable lattice models, when associated with rather strong 
symmetry requirements and with [9] or without [10] the so called "tetrahedron property" . 

However, this kind of consistency is very sensitive on the role of spectral parameters: the 
coefficients of the model have specific roles and are interdependent. This means, for example, 
that the normal set of linear transformations would completely destroy any such association. 
Unfortunately imposing this sort of interdependency adds enormous complications for a 
generic search program. 

1.3 Plan of the paper 

We start from lattice relations defined on an elementary square of the 2D lattice. We 
assume the relation is at most quadratic in the dependent variables. Due to the reasons 
discussed above we have chosen factorization as the first selecting criterion: we impose as 
a factorization requirement that (at least) a linear factor can be extracted after 2 steps of 
evolution. Solving the resulting equations produces a list of maps, as described in Section [3j 

There are natural symmetries between the models. For example two models may be related 
by a translation of the variables. Maps related by reflections or rotations can also be omitted, 
and this allows to reduce the list to 80 cases. 

We then analyze all 80 cases with an algebraic entropy calculation, as explained is Section 
HI This yields a classification of the models according to their degree growth. The vanishing 
of entropy (polynomial growth) is the integrability detector we use here. However, since 
our search condition is just the factorization of one linear factor after two iterations, many 
models with exponential growth are still included, which may contain integrable models, 
when further constraints are introduced. Our final list contains just the different "parents" 
without further analysis on their possible integrable "descendents" . 

In section \5\ we analyze in with some detail one particularly interesting multiparametric 
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Figure 1: The map is defined on an elementary square of the 2D lattice. 

model, and show that it also verifies the condition of consistency around the cube, thus 
proving its integrability. This model does not have the symmetries used in [9], and thus does 
not appear in the list given there, but as special cases it contains both the discrete modified 
and Schwarzian KdV. 



2 Factorization in 2D 
2.1 The map 

We consider maps defined on the cartesian 2D lattice by relating the four corner values of 
an elementary square (see Figure d]) with a multilinear relation: 

k XX[!]X[2]X[X2] + hxX[i]X[ 2 ] + hxX[i]X[i2] + h XX[2]X[ 12 ] + h X[l]X[ 2 ]X[ 12 ] 
+ p 1 XX {1] + p 2 X[!]X[ 2 ] + P3 X[2]X[t2] + Pi X[ 12 ]X + p 5 XX [2 ] + p 6 X[i]Z[i 2 ] 

+ nx + r 2 X[i] +r 3 X[2] +r 4 X[i2] + u = Q(x, X[i], X[ 2 ], X[i 2 y, «i, a 2 ) =0. (1) 

Here x n , m is the dependent variable at a corner and we have used a shorthand notation, in 
which only the shifts with respect to the base point at lower left is indicated in a subscript 
in square brackets: x nim = x 00 = x, x n+ i, m = x w = X[q, x n>m+ i = x QX = x [2 ], x n+ltm+ i = 
Xn = X[i 2 ]. (We will use indifferently these three notations). As mentioned above, in the 3D 
consistency approach an important role is played by the spectral parameters a s , associated to 
specific directions, and they appear in the coefficients k, Zj, Pi, r», u. In this paper, however, 
these spectral parameters are ignored (except in Section [5]), since we only work on the 2D 
lattice. 

Dynamics is defined by Q(x ntTn , x n+ i >m , x n>m+ ±, # n +i,m+i) = 0, using the multilinear Q of 
and this allows well-defined evolution from any staircase- like initial condition, up or down, 
since we can solve for any particular corner value in terms of the others, see Figure [2j 

The canonical examples of such maps are the following: 
• Lattice KdV 

(Pi ~P2 + x i -xi )(pi +P2 + x 00 - x u ) =p\ - pi, 
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or after translation x n>m = u n>m + p\ii + p 2 m 

Ooi -Uio)(«oo -Mn) =p\ -pi, (2) 

• Lattice MKdV 

Pi(x 00 xoi - x w xu) = P2(x 00 x w - x 01 xn), (3) 

• Lattice SKdV 

(x 00 - x w )(x i - xn)pl = (x 00 - x i)(x w - x X \)p\- (4) 
2.2 Keeping track of factors 

When a rational map is iterated some factors often get canceled "silently" in the resulting 
rational expression. This can be observed by comparing the degrees of the terms to the 
generic case. However, the best way to keep track of factors it to write the rational map 
as a polynomial map in a projective space. This was also the method that in the ID-case 
showed clearly the connection between singularity confinement and reduction in the growth 
of complexity [1]. In the ID-case it turned out that the amount of cancellation had to be 
such that the degrees only grow polynomially (oc n k ), while for nonintegrable systems we 
have exponential growth (p n with p > 1). The same happens for 2D-systems[T5l fT3] . 

Recall that the map we are considering is given by a multilinear expression 

This equation can be homogenized by substituting x n ^ m = v n>m / / n>m and taking the numer- 
ator. We assume that Q does not factor and that it depends on all the indicated variables. 
In the numerator of the homogenized Q we isolate i> n +i,m+i an d / n +i,m+i : 

^4 {vn,m ,Vn+l,m ^n,m+l i/n.m ■ ) fn-\-l,m,fn,m-\-l ) ^71+1,772+1 ~~ l~~ ( y n,m i^n+ 1 ,m ,^71,771+1 ifn.m i/n+l ,m ■,fn,rn-\-l ) ./n+l, 772+1 
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Figure 2: Fundamental evolutions on a square lattice. 
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Figure 3: Maximum degrees with staircase and corner initial states for quadratic relations. 



and then define the projective map as 



v n+l,m+l 
fn+l,m+l 




(5) 



The polynomials A, B are both of degree 3 in the indicated variables, and they cannot have 
common factors, since Q was assumed irreducible. 

For practical reasons we consider in this paper only quadratic maps, i.e., in (TjQ) we take 
k = 0, U — 0. In that case the projective map is given by 



where z = v or /, since they have the same degree. 

For the initial values we may take any staircase like configuration and in the rational repre- 
sentation take arbitrary x's at each point. In the projective representation we may also take 
arbitrary t>'s at all initial points, but we should use the same /, because projectivity only 
adds one free overall factor. In the quadratic case this means the cancellation of one extra 
term at each step, so that the sequence of maximal degrees (without any factorization) is 
1, 2, 4, 9, 21, . . . which corresponds to the asymptotically exponential growth (1 + v2)*. 

For the search part of our factorization study we will only consider the first few steps in 
the iteration of (JSJ). In principle two different initial configurations are often used, staircase 
and corner, the default degrees in these two cases are given in Fig. [3j The interesting 
factorization, that can be used to predict integrability, happens at the point where the 
default degree shown in Fig. [3] is 9 or 7, respectively. In the search part we use the corner 
configuration, because the total degree in then smallest. In the degree growth analysis we 
have used the staircase configuration. 




-[piVV m f [2 ] +P2V[ 1 ]V[ 2 ]f +P5VV[ 2 ]f[i] 

+n v /[i] f[ 2 ] + r 2 v {1] f [2] f + r 3 v [2 ] f[i] f + uf f {1] / [2] ] 

PZ V[2] f[l] f+P4V /[i] f[ 2 ] + p 6 U[i] f[2] f + r 4 / /[I] f[2] ■ 



(6) 



From this we see immediately that the default degree growth is: 



deg(z n+ i jm+ i) = deg(z n+ i tin ) + deg( ) + deg^m), 



(7) 
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2.3 Factorization of known models 



Consider the discrete KdV model given by 




n,m 




The homogenized Q is now 




n,m 



and the projective map 




(8) 



In the corner case we take f 0jTn = fn,o — foo an d after canceling one / o a ^ every step one 
find the first interesting GCD (greatest common divisor of the list of polynomials) at (2,2): 



Exactly the same GCD is found from the stair configuration, with initial values restricted 
by f n ,m = foo for n + m = 0,1. The remaining (different) factors of t>22,/22 are rather 
lengthy, in the corner case they are of degree 5 in the initial values and have 23 and 13 
terms, respectively, and for the staircase initial configuration the degree is 7 and the number 
of terms 112 and 76, respectively. 

The simplest linear model + xm + a; pi + x = the cancellations are so strong that the 
map stay linear, similarly for the quadratic model X\i2\X + axmxpi = the projective map 
stays quadratic. 

In [9] a list of lattice models having the "consistency around a cube" property was given. 
The divisors for the quadratic cases, using corner configuration, are as follows 

(HI) map: (x 00 - x u )(x 10 - x 01 ) + j3 - a = 0, 
factor (v w — Voi) 2 - 

(H2) map: (x 00 - x n )(x w - x m ) + {(3 - a)(x 00 + x w + x 01 + x n ) + (3 2 - a 2 = 0, 
factor: (v 10 - v i + (a - P)foo)(v 10 - v 01 - (a - 0)foo). 

(H3) map: a(x 00 x 10 + x 01 x n ) - (3(x 00 x i + x 10 x n ) + 5(a 2 - (3 2 ) = 0, 
factor: (v±oa — voifl)(voia — Vxo/3). 

(Al) map: a(x o + x 01 )(x 10 + x n ) - f3(x 00 + x w )(x i + x n ) - 5 2 a(3(a - (3) = 0, 
factor: (v 10 - v 01 + 5 (a - /3)/oo)(«io - v 01 - S(a - (3)f o)- 

(Ql) map: a(x 00 - x i)(x 10 - x xl ) - (3(x 00 - Xi )(x i - in) + 5 2 a(3(a - 0) = 0, 
factor: (v w - v 01 + 6 (a - /3)/oo)(«io - v o\ ~ 6 (a - (3) foo). 



GCD(w 2 2, /22) = (^oi - ^io) 2 - 
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(Q2) map: a(x 00 - x i)Oio - - Pi^oo ~ x 10 )(x 01 - x lt ) 

+ a (3 {a - P)(xqq + x 10 + x 01 + in) - a(3(a - (3){a 2 - a/3 + (3 2 ) = 0, 

factor: (v 10 - v 01 ) 2 - 2(a - /3) 2 /ooOio + «oi) + (a - /3) 4 /oo- 

(Q3) map: ((3 2 - a 2 )(x 00 Xn + x 10 x 01 ) + (3(a 2 - l)(x o£io + x 01 x u ) 

- a(p 2 - l)(x 00 x 01 + x w x n ) - S 2 (a 2 - f3 2 )(a 2 - l){[3 2 - l)/(4a/3) = 0, 

factor: 4af3(v Q ia - v m f3)(v m f3 - v w a) + 5 2 {a 2 - /3 2 ) 2 /oo- 

Note that for H1-H3, Al, Ql we get two linear factors while for Q2, Q3 the quadratic factor 
is irreducible (if 5 = then the divisor for Q3 does factor). Note also that the common 
factors can be used to study relationships between maps; for the above list we just observe 
that Al is obtained from Ql with x nm i— > (— l) n+m x nm . 

3 Search 

3.1 The method 

We have seen above that in all known integrable quadratic lattice maps the result at (2, 2) 
factorizes with a quadratic GCD (with the exception of linearizable models, that may have 
even more factors). In many models this quadratic common divisor factorizes, but not 
always. The search for a linear factor is computationally easier, and therefore this search is 
restricted to that, but hope to return to a search for model with irreducible quadratic factors 
later. The possibility of irreducible factors of degree higher than 2 is open, and no examples 
are known. 

In this search project we use the corner configuration, because computations are then simpler. 
Also, for computational simplicity, we use x rather than v, f, but to prevent accidental 
factorizations we proceed as follows: we calculate x±i, 221, £12 using a generic Q (which does 
not factorize), and substitute the obtained values into the equation Q(x±i, x±2, X21, X22) = 
and take its numerator (to be called Qn from now on). We then require that this Qn 
factorizes as Qn = (some polynomial in xoo, ^10, ^01, %2o, ^02, ^22) x [%oi + + soo]- Here 

we may assume that £01 has unit coefficient, because the other coefficients s can be rational 
in £00, furthermore, if £01 were zero but X\q not we could use n <-> m reflection. 

If we keep track of all the other variables except Xqo we have generically 



for some polynomial g, furthermore we know something about the exponents: a + ji + v + 
t + v < 7. Because of this observation we can make the following ansatz for the factorization 



1 



■1 




(9) 



a,fi,v=0 t,v=0 



s 



Q = PS: 



S =x i + sioxio + s 00 , (10) 

1 4 

P — x^x^x^x^qX^^Xoo, a, u, fi, t, v), (11) 

a,n,v=0 t,v=0 

where a + fi + v + t + v<6, and d, s are rational in xqq. 

Given the ansatz flTJ with the quadratic restriction k = 0, U = we compute the Qn (which 
has 15966 terms) and subtract the ansatz PS above. The equations are then formed by 
taking the coefficients of different powers of X22, £20, ^02, %io, ^01 with fixed ordering. Next we 
determined the 134 functions d in P by considering the leading terms of Q — P(xqi + . . . ) = 0. 
This is fully automatized using REDUCEjTT]. The remaining equations were then solved 
for the functions s and for the parameters of the map itself, pi,ri,u. The solution process 
branched a lot and required numerous separate computing sessions. 

4 Results 

If the computations led to a situation where the resulting map a) factorized, or b) did not 
depend on all corners, that branch was terminated immediately, but even then we got 125 
"raw" results. From this set we omitted maps that could be obtained by rotation or reflection 
from other maps (using translation by a constant if needed), this reduced the number from 
125 to 80. For all those remaining models we have calculated the algebraic entropy as 
explained in the following section. 

4.1 Growth patterns and algebraic entropy computations 

Suppose initial data is given on a line which allows the determination of the values at all 
points of the lattice, for example a regular diagonal staircase. The multilinear relation 
allows us to define evolutions in the following way: we iterate the relation by calculating the 
values on diagonals moving away from the initial staircase, as in the previous section. After 
cancellations we get a sequence of degrees d n . 

We may in this way define four fundamental evolutions, corresponding to initial data given 
on diagonals with slope +1 or —1, and evolutions towards the four corners of the lattice. We 
denote them "NE, SE, SW, NW" by the orientation of the evolutions (towards North-East, 
South-East, and so on), see Figure [2j To each evolution we associate an entropy with the 
definition inspired by the 1-dimensional case [SJ [T21 [TBI [T5] : 

e = lim i logK). (12) 

n— >oo 

These entropies always exist [12], because of the subadditivity property of the logarithm of 
the degree of composed maps. 
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A full calculation of iterates is usually beyond reach. We can however get explicit sequences 
of degrees as explained in [T3] . 



Suppose we start from initial values distributed on a diagonal, containing v vertices V\, . . . , V v . 
For each of these v vertices, we assign to the dependent variable x n ^ m an initial value of the 
form: 

a k + (3 k t 



k=l...v (13) 

where «o, fio and a k , f3 k , (k = l..u) are arbitrary constants, and t is some unknown. We then 
calculate the values of x at the vertices which are within the range of our set of initial date. 
These values are rational fractions of t, whose numerator and denominator are of the same 
degree in t, and that is the degree we are looking for. 

The next step is then to obtain "degree growth", i.e., to extract the value of the entropy 
from the first few terms of the sequence {d n }. One very fruitful method is to introduce the 
generating function of the sequence of degrees 
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(14) 



fc=0 



and try to fit it with a rational fraction. The remarkable fact is that this works surprisingly 
well, as it did for maps, although we know that it may not always be the case [H]. This 
means that we can often calculate the asymptotic behavior measured by (TT2]) from a finite 
beginning part of the sequence of degrees. 

We have done this calculation for all four evolutions of the 80 cases we got. The 4 x 80 
calculations resulted with a number of different sequences of degrees, which could all be 
fitted with rational generating functions. Some have linear growth, some are quadratic, 
some have exponential growth. The sequences are given in the following tables, containing 
the beginning of the sequences, the generating function, and the numerical value p of the 
growth p n 

Linear growth: 





Sequence of degrees 


Generating function 


Growth 


k 


1,2,3,4,5,6,7,8,9,10,... 




1 


h 


1,2,3,5,6,8,9,11,12,14,... 


1 + S + S A 


1 



Quadratic growth: 



(12 



Q3 



Sequence of degrees 



1,2,3,5,7,10,13,17,21,26, 



1,2,4,6,9,12,16,20,25,30,. 



1,2,3,5,7,11,14,20,24,32,. 



1,2,3,5,8,12,17,23,30,38,. 



1,2,4, 7,11,16,22,29,37,46, 



Generating function 

,'2,'i 



l-s z 



(s+l)(l-s) 3 

J_ £> £■ _L _1_ o' 



(s+l) 2 (l-s) 3 



Growth 
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Exponential growth: 





Sequence of degrees 


Generating function 


Growth 


ei 


1,2,3,5,8,13,21,34,55,89,... 


l+s 
1— s— s 2 


1.618 


e 2 


1,2,4,7,12,20,33,54,88,143,... 


(1_ S )(1_ S _ S 2) 


1.618 


e 3 


1,2,3,5,9,16,28,49,86,151,... 


1 

l-2s+s 2 -s 3 


1.755 


e 4 


1,2,3,5,9,17,32,60,112,209,... 


l-2s+s 3 -s 4 


1.867 


e 5 


1, 2, 4, 7, 13, 24, 45, 84, 157, 293, . . . 


1 

l-2s+s 3 -s 4 


1.867 


e 6 


1,2,4,8,15,28,52,97,181,338,... 


^1-f-SJ^l J 


1.867 


e 7 


1, 2, 3, 5, 9, 17, 33, 65, 129, 257, . . . 


1— s— s z 
fl-2»|(l-5) 


2. 


e 8 


1, 1, 3, 5, 11, 21, 43, 85, 171, 341, . . . 


(l+«)(l-25) 


2. 


e 9 


1, 1, 3, 6, 12, 24, 48, 96, 192, 384, . . . 


1— s+s 
l-2s 


2. 


eio 


1, 2, 4, 7, 14, 27, 54, 107, 214, 427, . . . 


(1- S )(1-2 S )(1+ S ) 


2. 


en 


1, 2, 4, o, lb, oz, o4, izo, zoo, olz, . . . 


l-2s 


Z. 




1 9/1 8 1fi ^9 1 QQ 97/I t^fifi 

1 , Z, 4, o, 10, OZ, 00, loo, Z I 4, 000, . . . 


(l +s )(l- s )' 
l-3s+* 2 +3s :i -2s 4 -s G 


Z.UO ( 


ei3 


1,2,4,8,16,33,68,141,292,605,... 


(l-s)( S +lJ 
l-2s-s 2 +2s 3 -s 5 


2.071 


ei4 


1, 2, 4, 8, 16, 33, 69, 145, 305, 642, . . . 


1-s-s 4 
(l-s)(l-2s-s 4 ) 


2.107 


eis 


1, 2, 4, 8, 17, 37, 82, 183, 410, 920, . . . 


(l-s)(l-2s-s 2 +s 3 ) 


2.247 


eie 


1, 2, 4, 9, 20, 45, 101, 227, 510, 1146, . . . 


(l-s)(s+l) 
l-2s-s 2 +s 3 


2.247 




1, 2, 4, 9, 21, 50, 120, 289, 697, 1682 . . . 


1-s-s 2 
(l-s)(l-2s-s 2 ) 


2.414 



In the previous table, the value e* is given for reference, it corresponds to the absence of 
factorization, and is an upper bound. The factorization condition we have used is a rather 
mild one, it says that the fourth number should be 8 or less. (The pattern ei6 does not have 
this cancellation, it arises in a model, which has cancellations in two direction, but not in 
the other two.) 

From these tables we also observe, that the sometimes a map with asymptotically exponential 
growth starts with as slow growth as is seen for linear growth, compare the patterns of 
ei, e 3 , e 4 , e 7 , e 8 with Z 2 . 

4.2 Classification into parents and descendents 

We are not going to give the full listing of the 80 models obtained, because the solution 
method is insensitive to subcase dependency. The criterion was only that there is one linear 
factor at position (2, 2) and this does not yet guarantee integrability. It is possible that one 
model is obtained from another one by specializing the values of the parameters, giving a 
notion of descendent. This process provides a partial ordering of the models, the descendents 
of a model having a smaller entropy. It may happen that a model with non- vanishing entropy, 
thus a priori non integrable, has a descendent with vanishing entropy. 

The list that is presented next is comprehensive in the sense that any integrable case with a 
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linear factor at (2, 2) will be a subcase of one of the presented models. The list is rather short 
and therefore further studies of subcases with more factorization will not be overwhelming, 
but not included here. 

The classification is up to transformations of the type 

x nm i-> x nm + a(n - no) + b(m - m ) + c, 

where the a, b terms can only be used if the result is n, m independent. We would also like 
to note that innocuous reparametrizations can change the factorization properties, because 
computer factorization is normally over integers. Examples of this can be seen in Cases 2 
and 5. 

The ancestral models are as follows: 
Case 1 

This is a homogeneous model, where one quadratic term is missing, all others have arbitrary 
coefficients: 

^oo^ioPi + x 00 x u P4 + x w x 01 p 2 + x 10 x u p 6 + x 01 x n p 3 = (15) 

The degree sequences in the four directions are eis, exe, eie, ei5. The missing term above is 
Xoo^oi) but it could be rotated into any other side of the elementary square. Since the model 
is not rotationally symmetric it can have different growth patterns in different directions. 
The pattern e 1& corresponds to growth without any cancellations at (2,2), but e 15 , which 
has some cancellations, is obtained in the direction of the test. As for its subcases, we 
just mention that if p 3 = 0, i.e., the opposite side xio^n is also missing, then the map is 
linearizable. 

Case 2 

This is also homogeneous, but it has "symmetric cross" xoo^ii + a^io^oi and all other coeffi- 
cients are arbitrary 

xooXioPi+xooXo 1 p 5 (p 1 p 3 + p 2 ) + (x 00 xn+x 1 QXo 1 )p2 + x w xnp fi + xo 1 xnP3(P5P(i-P2) = (16) 

This has quadratic growth q§ is all directions. This result illustrated the role of parametriza- 
tion in factorization: there are several ways to reparametrize P3,Ps so that combinations 
P5(PiP3 + P2) and Pzip^Ps — P2) look simpler, but this could easily lead to one quadratic 
rather than two linear factors at the corner. 



The next four models have a free nonhomogeneous constant term. 
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Case 3 



Here there are completely arbitrary linear and constant terms, while in the quadratic part 
two opposite sides are missing (a^o^io and XqiXu) 

x 00 x 01 p 5 + XooXnPi + x w x 01 p 2 + x 10 x u p 6 + x 00 ri + x w r 2 + x 01 r 3 + x n r 4 + u = (17) 

The growth pattern is en in all directions, so as it stands the model is not integrable. Among 
its subcases we would like to mention 

(P3X Q1 + x o)(P3Xh + xio) + (P3X01 + x w )r 3 = (18) 

and 

(x 00 - x i)(x 10 - x u ) + (x 00 - x u )r 4 + (x i - x w )r 3 + u = (19) 
both integrable, with growth patterns q 3 in all directions. 

Case 4 

Here is the homogeneous part the cross is missing, and other terms have arbitrary coefficients, 
except that there is a relation between the coefficients: 

^llXioPe + X U X ip 3 + XioXooPl + 34)1^00P5 ^ 2Q ^ 

+ Xup 3 p & ri + x w p 6 r 2 + x 01 p 3 r 3 + x 00 (-piP 5 ^4 + Pi^3 + ^5^2) + u = 0. 

The parameters have been scaled to simplify the model, but still nicer forms of writing it 
may exist. The growth pattern is again en in all directions. This contains as subcases the 
following integrable ones: 

(x ox 01 + xiox n )p 5 + (x 00 x w + x 01 x u )p 3 + u = 0, (21) 

(x 00 - x u )(x 01 - x w ) + ri(x 00 + zio + Xoi + zn) + u = 0, (22) 
which include H1,H2,H3 of [9] 

Case 5 

The following model is integrable, with pattern q 5 is all directions: 

(x n + x 00 )(a;io + x i)p 2 + (x n + a; i)(x o + x w )p 3 + p 3 p 2 {p 3 + p 2 )u 2 = 0, (23) 

or in another way of writing 

(x 11 x w +XoiXoo)(p 2 +p 3 ) + (xi 1 x 01 +x 10 x 00 )p 3 + (x 11 Xoo+x w x i)p 2 +p 3 p 2 (p 3 +p 2 )u 2 = 0. (24) 

The above parameterization of the constant term is necessary in order to have two linear 
factors at (2,2), otherwise we would get one quadratic factor. This case becomes Ql with 
scaling x nm 1— > (— l) n ~ n °x nm (it is not separately listed in j9], but the situation is similar to 
Al). 
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Case 6 



XuXoo + XioXoi + (Xn^oi + Xi X 00 )p3 - (xnxio + x ix 00 ){p 3 + 1) ^ 
+ (x u - x 00 )r 4 + (xio - x i)r 2 - (u(p 3 + 1) + r 4 )(up 3 + r 4 ) + ur 2 = 0, 

which is integrable with growth q§ in all directions. If r 4 = r 2 = we get Ql of [9]. In fact, 
if p 3 7^ 0, — 1 we can apply an n, m dependent translation to put r 2 = r 4 = 0. 



The remaining maps have linear and constant terms, but the constant term depends on 
the linear terms. There are also interesting relationships between the parameters of the 
homogeneous terms. There may be other interesting representations obtainable with the 
transformation Xy — > [ix^ + v. 



Case 7 



x\\x 1{i v\ + xoiXooPe + x X \Xq\PzV\ + xvdX^p 3 p x 1 pI 



+ (x n P3 + x 10 p 3 + x ipa + XooPeK + r t 
with growth en in all directions. 



( 26 ) 



Case 8 

xnx w pl + x 01 x 00 pl + XnSoiPi ^(Pa - 1) + ^lo^ooPiPefe - 1) + (zn^oo + ^io^oi)Pe 
+ (xuPe + x w pe + xoipz + x 00 p 3 )r 4 + r\ = 
also with growth en in all directions. 



Case 9 

(xuXqo + x 10 x i) + psixuxoi + x w x 00 ) + pe^n^io + x 01 x 00 Pe 1 (p3 - I) 2 

+ ^3(^6 -P3 + l)(a?oi +^00) + Pe(pe + l)r| = 0. 
This has growth e 44 in all directions. In the subcase p^ = p 3 — 1 it becomes the integrable 
case Al of [5], on the other hand, if p§ 7^ p 3 — 1 it can be translated into the more symmetric 
form 

(xnXoo + XiqXqx) + P^XuXqx + X10X00) + PqXux w + x 01 XooP6 1 (P3 - I) 2 , s 

2 ' J 

+ r 4 (xn + xio + £01 + ^oo) + r 4 = 0. 



Case 10 

OoO^lO + XqiXu) + (xooXn + X lX W )p 3 + (j9 3 - 1)(X 00 X 1 + XioXu) ^ 

(x 00 - x Q1 + X10 - xn)r 4 + 7-4 = 0, 
also with growth e u in all directions. 



14 



Case 11 



(x 00 xu + X01X10 - x ix 00 )pa + Pix 00 x 10 + (xooPi + xio + xiip 4 )r 4 + r\ = (31) 
with growths e$, e±e, e^, in the for directions. 

Case 12 

^oo^ioPi + x 01 x u p 3 + ^oo^oi (P1P3 - 1) + ^oo^n + ^io^oi ^ 
+ (xoopx + xio + x 01 p 3 + xn)r 4 + r\ = 

with growths es, e±Q, e\Q, e\§. 

The known models referred to in section (12. 3p . that have linear factors, all appear in our 
analysis. 



5 Some integrable cases 



Due to the search condition (one linear factor at (2,2)) the results are not comprehensive 
as far as integrable cases are concerned. We can say that they must all be subcases of the 
cases enumerated above. 



Case 2 is integrable as it stands. We can also write it as 

^oo^ioCi + ^oo^oiCs + (^oo^ii + ^io^oi)^ + + %m%iiC3 — 



(33) 



and the statement is that this lattice map is integrable for all values of the five parameters 
Cj. If we next consider integrability in the sense of consistency, then the parameters must 
have specific forms in terms of the spectral parameters associated with the lattice directions. 
The result is 

(xx [ij] + x [j] x )(ai-aj) + xx w (/% + af) / jj - xx (A + a^/^-x H x ^ ji + x^x jj = (34) 
where the parameters associated with the third direction are restricted by 



1 ai Pi 
1 a 2 fa 
1 «3 Aj 



0. 



(35) 



An immediate way to resolve this constraint is to take 

fip} + v, f3 i = pp 2 i +a. 



ft; 



The various homogeneous KdV equations are obtained as subcases (with one spectral pa- 
rameter), the modified KdV of ([3]) is obtained with Ji — Pi, p — — 1, fi — u — a — and the 
Schwarzian KdV of (TjJ with ^ = —pj, /i = —1, p = v = a = 0. Furthermore the generalized 
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form of [IB] is included with 7, = — (pi + a)(pi + b), pL — — 1, v — 0, p — — (a 2 + b 2 ), a = a 2 b 2 . 
By suitable scaling of type x nm 1— > A n ~ n °B m ~ m ° x nm one can reduce this model to Q3 with 

5 = 0. 

It is perhaps worth mentioning that if some of the q vanish we may get asymmetric models 
with different growth in different directions, for example 

a^oo^ioci + xoo^oiCs + (x ox n + xi x i)c 2 = 0, (36) 

has growth patterns g 2 , Q5 , qi , Q5 , and 

^oo^ioci + x 10 x u c 6 + x 01 x n c 3 = (37) 

grows as (fc, <?4, ^4, <?5- It is not clear how such special cases carry over to the consistency 
approach. 

As an interesting model worth further study we would like to present ( fT9l 

(x 00 - £oi)0eio - X11) + (^00 - xu)r A + (x i - x w )r 3 + u = 0. 

It is integrable with growth patters ^3 in all directions. But this model is not symmetric 
and therefore its consistency formulations is problematic: How should the maps on the other 
sides of the cube be oriented, or do we perhaps need entirely different maps there? 

6 Conclusions 

We have analyzed the factorization process for a class of two-dimensional lattice models 
("quad models"). By solving the resulting equations we obtained a list of models. We then 
conducted an entropy analysis on the results. This provided an ordering among the models, 
by increasing complexity. 

This list of minimally factoring models contains both models with vanishing entropy, and 
models with different non-vanishing entropy. These models with non-vanishing entropy may 
contain integrable models as subcases, but this has not been comprehensively analyzed here. 

More can be done along the lines we have followed: one could insist on stronger factorization 
requirements, or one could start from the most general quartic defining relation, rather than 
restricting on quadratic ones as we did here. Both of these is beyond the scope of this paper. 
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